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Abstract 



Theories unifying gravity and other interactions suggest the possibility of spa- 
tial and temporal variation of physical "constants" in the Universe. Detection 
of high-redshift absorption systems intersecting the sight lines towards distant 
quasars provide a powerful tool for measuring these variations. We have previ- 
ously demonstrated that high sensitivity to the variation of the fine structure 
constant a can be obtained by comparing spectra of heavy and light atoms (or 
molecules). Here we describe new calculations for a range of atoms and ions, 
most of which are commonly detected in quasar spectra: Fe II, Mg II, Mg I, 
C II, C IV, N V, O I, Al III, Si II, Si IV, Ca I, Ca II, Cr II, Mn II, Zn II, Ge II 
(see the results in Table 3). The combination of Fe II and Mg II, for which 
accurate laboratory frequencies exist, have already been used to constrain a 
variations. To use other atoms and ions, accurate laboratory values of frequen- 
cies of the strong El-transitions from the ground states are required. We wish 
to draw the attention of atomic experimentalists to this important problem. 

We also discuss a mechanism which can lead to a greatly enhanced sensitiv- 
ity for placing constraints on variation on fundamental constants. Calculations 
have been performed for Hg II, Yb II, Ca I and Sr II where there are optical 
transitions with the very small natural widths, and for hyperfine transition in 
Cs I and Hg II. 
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I. INTRODUCTION 



Possible variations of the fundamental physical constants in the expanding Universe are 
currently of particular interest because of the implications from unified theories, such as 
string theory and M-theory, that additional compact dimensions of space may exist. The 
"constants" seen in our three-dimensional subspace of the theory will vary at the same rate 
as any change occurring in the scale lengths of the extra compact dimensions (see, e.g. JH-Q). 
Gas clouds which intersect the sightlines towards distant quasars produce absorption lines. 
These absorption systems present ideal laboratories in which to search for any temporal or 
spatial variation of fundamental constants by comparing the observed atomic spectra from 
the distant objects with laboratory spectra (see, e.g. || and references therein). 

The energy scale of atomic spectra is given by the atomic unit I p-. In the non-relativistic 
limit, all atomic spectra are proportional to this constant and analyses of quasar spectra 
cannot detect any change of the fundamental constants. Indeed, any change in the atomic 
unit will be absorbed in the determination of the red shift parameter z (1 + z = ^7, u' is the 
red-shifted frequency of the atomic transition and to is the laboratory value). However, any 
change of the fundamental constants can be found by measuring the relative size of relativistic 
corrections, which are proportional to a 2 , where a = e 2 /he is the fine structure constant ||. 

It is natural to search for any changes in a using measurements of the spin-orbit splitting 
within a specific fine-structure multiplet, and indeed this method has been applied to quasar 
spectra by several groups. However, whilst this method is appealing through its simplicity, it is 
possible to improve on its efficiency Also, it may even give incorrect results, since corrections 
of higher order than a 2 are ignored. An order of magnitude sensitivity gain can be achieved by 
comparing transition frequencies of heavy and light atoms (or molecules) In this paper 

we extend that previous work, presenting the results of calculations of the dependence of the 
transition frequencies on a for many atoms and ions where data exist for quasar spectra (see 
Table 

Other possibilities for measuring changes in a include comparisons of different optical tran- 
sitions, such as s — p and p — d, in the same atom or molecule, or comparisons of microwave 
transitions in molecules which contain rotational and hyperfine intervals. We have also calcu- 
lated the dependence on a of some atomic microwave and optical frequency standards which 
could be used for laboratory searches for a variations. For example, one can compare the 
Hg II E2 transition uj = 35514 cm -1 with any narrow line of another atomic or molecular 
transition having approximately the same frequency. Small frequency differences (which do 
not require very precise absolute calibration) can be measured very accurately. The mercury 
frequency has a very large relativistic shift which also has a negative sign (usually the shift is 
positive). Therefore, this frequency difference has a very strong a-dependence. 

Finally, there is an interesting possibility for studying transitions between "accidentally" 
degenerate levels in the same atom or molecule. There are several practically degenerate 
levels of different electron configurations in rare-earth and actinide atoms. Of course, there 
are many more possibilities in molecules where there are vibrational and rotational structures. 
The relativistic corrections to the different energy levels are different and can exceed the very 
small frequency corresponding to the transition between "degenerate" states by many orders 
of magnitude, i.e. a tiny variation of a can change the frequency significantly. Also, there 
is an interesting dependence on the nucleon mass if the rotational, vibrational and hyperfine 
structures are involved. In this case one can measure time dependence of the nucleon mass 
which is a function of the strong interaction constant. The main problem here is to find a 
narrow-width transition. 
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Note that we present all results in this paper assuming that the atomic unit of energy ^p- 
is constant. 



II. THEORY 



A. Semi-empirical estimations 

Let us start our calculations using simple analytical estimates of the relativistic effects in 
transition frequencies. First consider the relativistic corrections to the frequency of an atomic 
transition in a hydrogen-like atom. The relativistic correction to the energy level is given by 
(see, e.g. §) 

me*Z*(Za)\ 1 _3_ 
n ~ 2h 2 S' + l/2 4n J ' [) 

where Z is the nuclear charge, n is the principal quantum number and j is the total electron 
angular momentum. This value of the relativistic correction can be obtained as an expectation 
value (V) of the relativistic perturbation V, which is only large in the vicinity of the nucleus. 
Therefore, the relativistic correction A is proportional to the electron density near the nucleus 
| \& (r < -|) | 2 oc ^3 (a is the Bohr radius, % is the size of the hydrogen- like ion) . For an external 
electron in a many-electron atom or ion the electron density near the nucleus is given by the 
formula (see, e.g. |J) obtained in the semi-classical approximation {n 3> 1) 

l*(r<|)|vg, (2) 

where Z a is the charge "seen" by the external electron outside the atom, i.e. Z a — 1 for 
neutral atoms, Z a = 2 for singly charged ions, etc.; v is the effective principal quantum 
number, defined by E n = — ffr-^f, where E n is the energy of the electron. For hydrogen-like 
ions v = n, Z a = Z. Thus, to find the single-particle relativistic correction, we should multiply 
A in Eq. ([I]) by the ratio of |\l/(r < -|)| 2 in the multi-electron ion and hydrogen-like ion. The 
result is 



me'Zl (Za f 
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The second term in the square brackets is presented to provide a continuous transition from the 
hydrogen-like ion Eq. (^) to the multi-electron ion Eq. (^). In multi-electron ions (Z ^> Z a ) 
this term is, in fact, a rough estimate based on the direct calculation of (V). We should 
neglect this small term since there are more important many-body corrections. 

We see that the relativistic correction is largest for the Si/ 2 and pi/ 2 states, where j = 1/2. 
The fine structure splitting is given by 

A is = E(p 3/2 ) - E( Pl/2 ) ~ -A(p 1/2 )/2 ~ -A(p3 /2 ). (4) 

In quasar absorption spectra, transitions from the ground state have been observed. There- 
fore, it is important to understand how the frequencies of these transitions are affected by 
relativistic effects. The fine splitting in excited states is smaller then the relativistic correction 
in the ground state, since the density of the excited electron near the nucleus is smaller. As a 
result, the fine splitting of the Sl-transition from the ground state (e.g., s—p) is substantially 
smaller than the absolute shift of the frequency of the s—p transition. The mean energy of 
the p-electron is defined as 

E(p) = ^E(p 3/2 ) + l -E(p 1/2 ) ~ E n (p) - ^A ls , (5) 

where E n is the non-relativistic energy. Therefore, the relativistic shift of the mean s—p 
transition frequency is given by 



A(p-s) ~ --A Zs - A(s 1/2 ). 



(6) 



The formulae (|3]) - (|6]) do not take into account many-body effects. For example, relativistic 
corrections change the self-consistent atomic potential. The many-body calculations discussed 
below show that the relativistic energy shift in atoms with one external electron can be 
approximately described by the following equation 



^(Za) 2 
v 



' C(Z,j,l) 



3 + 1/2 



(7) 



where C(Z,j, I) is different for different atoms and partial waves but does not depend on the 
principal quantum number. In many cases C(Z,j,l) ~ 0.6, although noticeable deviations 
from this value are possible. 

It is easy to explain the sign of the many-body effect. The relativistic single-particle 
correction increases the attraction of an electron to the nucleus and makes the radius of the 
electron cloud smaller. As a result, the direct Hartree-Fock potential, which is the nuclear 
potential screened by core electrons, becomes smaller at short distances. This decreases the 
binding energy of a valence electron. Therefore, many-body effects have an opposite sign to 
the direct effect. 

It is also easy to see why C(Z,j, I) does not depend on the energy of the valence electron. 
The effect of a potential change on short distances on the binding energy of an external 
electron is proportional to the density of this electron in the vicinity of the nucleus which in 
turn is proportional to 1/z/ 3 (see eq. (0)). On the other hand, the direct relativistic effect 
is also proportional to 1/z/ 3 for the same reason (see eqs. (@^J)). Therefore the ratio of two 
effects which is C(Z,j,l)/(2j + 1) is practically independent of the energy of the external 
electron. 

The accurate value of the relativistic shift can only be obtained from many-body calcu- 
lations. However, one can assume that C(Z,j,l) = 0.6 and use (0) for rough estimates of 
relativistic corrections. This usually gives better results than a single-particle estimate in eq. 
(0). For example, many-body calculations show that, as a rule, the relativistic correction to 
the energy is negative for si/ 2 and pi/ 2 states and positive for other states. This behavior is 
reproduced by (0) but not by (|3|). Apart from that, formula (|j) suggests that the correction is 
largest for s\/ 2 and p\/ 2 states and rapidly decreases with j while many-body calculations show 
that the correction for (i-states is sometimes bigger than for p-states. This is again reproduced 
by (|7|) where there is a strong cancelation between two terms in the case of p-states. Note that 
this complex behavior of the relativistic effects cannot be explained in terms of single-electron 
density at the origin and should be attributed to many-body effects. 

Let us see now how eqs. (|])-(0) wm t> e modified if the many-body correction C(Z,j, I) is 
included. Assuming the same value of C(Z,j,l) for pi/ 2 and p 3 / 2 states we have for the fine 
splitting 

A, s = E(p 3/2 ) - E(p 1/2 ) ~ 4 E %(j+T/2) ' (8) 

For the mean energy of the p-electron instead of (|^) we have 

E(p) = E n (p) + (2C - ^)A ls . (9) 

With C = 0.6 the second term in @ is numerically small. Hence, relativistic corrections move 
the pi/2 and ps/ 2 states in opposite directions, leaving the mean energy almost unchanged. 
Therefore, the relativistic shift of the mean s—p transition is mostly given by the energy shift 
of the s-state 



A(p-s) ~ -A(s 1/2 ). 



(10) 



The relative sizes of relativistic corrections are proportional to Z 2 , so they are small in light 
atoms. Therefore, we can constrain changes in a by comparing transition frequencies in 
heavy and light atoms. We stress that the most accurate and effective procedure includes 
all relativistic corrections and the analysis of all available lines (rather then the fine splitting 
within one multiplet only). We have not discussed here contribution of the Breit relativistic 
correction to the electron-electron interaction. It is not enhanced by the factor Z 2 and is 
much smaller then the contribution from spin-orbit interaction. Our numerical calculations 
have demonstrated that the contribution of the Breit interaction to the frequency of s-p 
transition A(p — s) is indeed negligible. 



B. Relativistic many- body calculations 

Accurate calculations of relativistic effects in atoms have been done using many-body 
theory which includes electron-electron correlations. We used a correlation-potential method 
Il0| for atoms with one external electron above closed shells and a combined configuration 



interaction and many-body perturbation theory method |TT] for atoms with several valence 
electrons. A relativistic Hartree-Fock (RHF) Hamiltonian was used in both cases to generate a 
complete set of single-electron orbitals. The value of the relativistic corrections were obtained 
by repeating the calculations for different values of a in the RHF Hamiltonian. 

The form of the the single-electron wave function we use explicitly includes a dependence 
on the fine structure constant a 

Mr) 7 = - ( ^ r )™^ r / r W ^ (11) 
r {iag(r) n n(v/r) jlm ) " {iL) 

This leads to the following form of the RHF equations 

f n {r) + -fn(r) - [2 + a\e n - V)]g n {r) = (12) 
r 

9n(r) ~ -9n(r) + (e n - V)f n (r) = 0, 
r 

where k = (— l)'+J +1 / 2 (j + 1/2) and V is Hartree-Fock potential: 

Vf = V d (r)f(r) - [ V exch (r, r')f(r')dr'. (13) 



The non-relativistic limit can be achieved by reducing the value of a to a = 0. 

For atoms or ions with one external electron above closed shells, the calculations begin 
with the V N ~ X approximation. A Hartree-Fock procedure is carried out initially for a closed 
shell ion with the external electron removed. When convergence for the core is achieved, the 
states of the external electron are calculated in the field of the frozen core. The main reason 
for this approach is the simplicity of the perturbation theory for calculation of the correlation 
corrections. It is well known that correlations are very important in many-electron atoms and 
should be included into the calculations to obtain accurate results. We do this by means of 
many-body perturbation theory (MBPT) and a correlation potential method [|ITJ . 

The correlation potential £ is defined as an operator which gives a correlation correction 
to the energy (ionization potential) of the valence electron 

Se a = (a\t\a), (14) 

The expectation value here is taken over the single-particle wave function of the external 
(valence) electron. Thus, £ is another non-local operator which can be included in (^) by 
redefining the non-local potential 



Vexchfar') -> V exch (r,r') + £(r,r'). 



(15) 



Single-electron states of a valence electron calculated in a non-local potential flTq ) are often 
called Brueckner orbitals. 

MBPT is used to calculate S. Perturbation expansion starts from the second order in 
residual Coulomb interaction. Dominating higher-order correlations can also be included by 



using the technique developed in Refs. [12]. However we found that for the few electron ions 
such as C IV and Si IV, second-order results are already of very high accuracy (see Table 
|). Therefore, we decided that the accurate calculations of higher order corrections are not 
needed. 

For some other atoms we introduce fitting parameters f v into the expression for the non- 
local potential in (0) to simulate the effect of higher-order correlations: V non -iocai = V exc h + 
/„£, where v = s,p or d. The values for f 8 ,f p and fd are chosen to fit the experimental 
data for laboratory value a = 1/137.036. In all cases the values of /„ are close to unity. The 
same f v have been used for calculations with varying a. This procedure works well because 
the accuracy of the results in the second order is already good and only a small correction is 
introduced by the fitting parameters. 

For atoms with more then one external electron we use the combination of the configuration 
interaction method with many-body perturbation theory |TT|. 



1. As for single-electron-above-closed-shells atoms, we start calculations from the RHF 
method in V 1 ^" 1 approximation. However, this starting approximation does not corre- 
spond to a closed-shell system and needs to be further specified. We do this in a very 
simple way. The contribution of an open shell to the Hartree-Fock potential is calcu- 
lated using the complete shell potential multiplied by the "occupation" factor nj (2j+ 1), 
where n is the actual number of electrons on that shell and j is the total single-electron 
momentum. Single-electron basis states are calculated in this Hartree-Fock potential 
with one external electron removed (V^ -1 potential). 

2. All basis states are divided into core states and valence states. Core states are frozen 
and included into calculations only via the effective potential of the core. Valence states 
are used as a basis for the configuration interaction method. Note that the definition of 
the core at this stage does not necessarily coincide with the core in RHF calculations. 
For example, for Fe II 3d^/2 and 3d 5 / 2 states are core states in the RHF procedure. 
But these states are also included into the configuration interaction, so they are valence 
states in the CI calculations. 

3. The effective Hamiltonian of the CI method is constructed. To include correlations 
between core and valence electrons we modify the effective Hamiltonian of the standard 
CI method by adding an extra operator £ 

if CI fjCI , y> 

n eff ~^ U eff + Z - 

The operator £ consists of two parts. £i is the one-electron operator which describes 
the correlation interaction between a valence electron and the core. £i is very similar to 
the correlation potential £ in ( jl5|) which we used for atoms with one external electron 
above closed shells. £ 2 is a two-particle operator which describes the effects of screening 
of the Coulomb interaction between valence electrons by the core electrons. 



It has been demonstrated in Refs. [£L1[] that the core-valence correlations are very im- 
portant and usually dominate over correlations between valence electrons. Thus, it is 
more important to include £ than to achieve the completeness of the basis in the CI 
calculations. 



4. The standard CI technique is used to diagonalize the matrix of the effective CI Hamil- 
tonian and obtain many-electron energies and wave functions. 

We need to apply this method to the atoms with many electrons in open shell. This can 
make configuration space and computation time very large. However, we do not need very 
high accuracy in the present calculations. Therefore, we have made a few simplifications to 
the method compared to Ref. [fil l. The effect of screening (E 2 ) is usually much smaller than 



one-electron correlations with the core (X^) and we neglect it. We also neglect subtracting 



diagrams in Y,i (see Ref. [[I ij) because there must be strong cancelations between screening and 
subtracting diagrams. Finally, we use a relatively small basis set with 3 to 4 single-electron 
basis states of each symmetry. This usually leads to a few hundred configurations. This small 
basis is not complete to high precision. However the results are reasonably good because 
correlations with the core are included (see Table [Tl|). To simulate the effect of incompleteness 
of the basis and of the omitted diagrams we introduce fitting parameters for S, similar to 
what we did for atoms with one external electron. The results for Cr II presented in Table [Tl] 
illustrate the effects of core-valence correlations, correlations between valence electrons and 
the effect of fitting. Note that only one fitting parameter was used to fit all energy levels. 

There are two contributions to the relativistic energy shift. The first is the direct rel- 
ativistic correction to the energy of a valence electron in the Hartree-Fock-Dirac equations 
(|12|). This correction can be found be varying a in eq. ( |T2D with fixed potential. There is 
also an indirect relativistic correction which appears because of the change of the core po- 
tential (including correlation potential) due to the relativistic effects in the core. Neither 
the correlation potential E, nor the standard Hartree-Fock potential depend on a explicitly. 
This dependence appears via the basis set of single-electron wave functions used to calculate 
both potentials since these wave functions had been obtained by solving Dirac-like equations 
(|Ti2|). Full-scale calculations, repeated for different values of a, were necessary to reveal this 
implicit relativistic behavior. The value of the indirect effect is not small and in some cases 
exceeds several times the value of the direct one. This indirect relativistic effect is essentially 
a many-body effect. It is responsible for the failure of the single-particle formula (|3|) which is 
unable to reproduce the value of the relativistic correction accurately. 



III. RESULTS AND DISCUSSION 

To find the dependence of frequencies on a we use the following formula for the energy 
levels within one fine-structure multiplet: 

E = E + Q!((-) 2 - 1) + Q 2 ((-) 4 - 1) + ^(LS)(-) 2 + if 2 (LS) 2 (-) 4 . (16) 

Here E , Ql and Q 2 describe the position of the configuration center, K\ and K 2 describe 
the level splitting within one configuration, L is the total orbital angular momentum, S is 
the total electron spin and ai is the laboratory value of a. We introduce an (LS) 2 term to 
describe deviations from the Lande interval rule. There are two sources of the (LS) 2 term: the 
second order in the spin-orbit interaction (~ (Za) 4 ) and the first order in the Breit interaction 
(~ a 2 = 5.3 x irr 5 ). 

The second-order spin-orbit interaction is larger for heavy atoms where we actually need to 
introduce the (LS) 2 term. Therefore, we first fitted the experimental fine structure intervals 
to find K\ and K 2 (numerical calculations give close values of K\ and K-i). Then we used 



numerical calculations for a = y7/8oti and a = y3/4ai to find the dependence of the config- 
uration center on a (coefficients Qi and Q 2 )- It is convenient to represent the final results in 
the form 



uj = uj + q 1 x + q 2 y, 



(17) 



where q x = Q x + #i(LS),g 2 = Q 2 + K 2 (LS) 2 , x = (ff - l,y = (^) 4 - 1 and u = 
Eq + i^i(LS) + ^(LS) 2 is an experimental energy of a particular state of the fine struc- 
ture multiplet. The parameters u>o, q x and q 2 for El transitions for many atoms and ions of 
astrophysical interest are presented in Table |TTJ. 

One can use these data to fit absorption systems in quasar spectra in order to measure or 
place upper limits on any variation of a. The maximum theoretical sensitivity comes from 
comparing the spectra of Fe II and Cr II since relativistic effects in both ions are large and 
have opposite sign. The effect here is about 20 times larger than the fine splitting for each of 
the ions. 

An analysis of the theoretical data reveals some interesting tendencies in the behavior of 
the relativistic corrections apart from being proportional to Z 2 : 

• Within a series of one-electron orbitals of a given symmetry, the relativistic energy shift 
is largest for the lowest orbital and decreases for orbitals of higher energies. This trend 
is supported by the semi-empirical consideration presented in section 1. The higher 
electron density in the vicinity of the nucleus gives larger relativistic effects. This also 
explains why the relativistic shift of the ground state energy is bigger when ionization 
potential is larger. 

• When a is changing towards its non-relativistic limit a = 0, one-electron energies of s 
and pi/2 states move up while energies of ps/ 2 and d states move down. The relativistic 
shifts of s and d states tend to be large while the energy shifts of p states are relatively 
small. Note that the single-particle consideration suggests that the relativistic shifts are 
large for s and pi/ 2 states and all energies move up when a is decreasing (see formulae 
(0-0). Accurate calculations give different behaviors due to indirect relativistic effects: 
relativistic corrections to the core orbitals change the electronic potential which in turn 
shifts the energy of the external electron. This effect is neglected in a naive single-particle 
formulae. The direct effect dominates over the core change effect for s-states. For p- 
states these two effects are close in magnitude but opposite in sign. The core change 
effect dominates for d-states. Thus, introduction of C(Z,j,l) in eq. (||) effectively takes 
into account the core change effect. 



These tendencies are illustrated by the results presented in Table [III]. For example, the 
relativistic shifts of frequencies of El transitions are negative in Cr II and Ge II and positive 
for other atoms. This is because these transitions correspond to a d—p one-electron transition 
in the case of Cr II and to a p — s transition for Ge II while for the other atoms the transitions 
are of the s — p type. The relativistic correction for Ge II is relatively small because the 
corresponding transition may be described as the transition from the ground p orbital to 
excited s orbital and because the relativistic corrections are small for both p states and excited 
states. Another example is the sharp increase of the relativistic effect from Ca II to Cr II. 
This is due to the coherent effect of two factors: bigger Z and bigger ionization potential for 
Cr II as compared to Ca II. 

As is apparent from the analysis above one should expect the biggest relativistic shift for 
the s — d (or d — s) transitions in heavy atoms. These transitions are not observed in quasar 
absorption systems but may be suitable for laboratory experiments. The natural line width for 
these transitions is very small. Futhermore, very precise measurements of the frequencies of 
many such transitions already exist since they are used as atomic optical frequency standards. 
We present in Table [TV] relativistic shifts of frequencies of some atomic transitions which are 
used or proposed as optical frequency standards. These include the strongly forbidden El 
transition in Ca I [13] and E2 transitions in Sr II [H], Yb II |nj and Hg II [IB"]. There are 



many other atoms which are being studied as possible frequency standards but which are not 
included in the table. These include Mg I PH , In II PH, Xe I PH, Ag I [g0|, etc. Note that 



the biggest relativistic effect is in Hg II. This is because of the d — s transition and a high 
value of Z. This makes Hg II the most interesting candidate for a laboratory search for a 
variations. 

The following limit of a variation was found in Ref. J7|: 

®- = +2.6 ±5.2 x 10- 16 (18) 
Assuming at/ at = 10 -15 one can get for the 5d l0 6s 2 S\/2 — 5d 9 6s 2 2 D 5 / 2 transition in Hg II 

Cj = 3Hz yr _1 , 



which should be compared with the natural linewidth limit of 1.8 Hz [16( . 

There are also ongoing laboratory searches for variations of a using microwave atomic fre- 
quency standards (atomic clocks) (see, e.g. ||22|| ). There are a number of microwave frequency 
standards which use ground state hyperfine structure (hfs) intervals of atoms or ions (see, e.g. 
review J23|). These include hfs of Rb, Cd + , Cs, Ba + , Yb + , Hg + f23fl . Here again the biggest 



relativistic effects are in the Hg + ion. It is convenient to present the a-dependence of hfs 
constants in a form similar to (fF, 



A hfs = (-) 2 (A + qx), (19) 
ai 

where Aq is the hfs constant for a = a.\. Its value was recently measured to very high precision: 

A (Hg + ) = 40507347996.84159(14)(41)Hz [fH]]. 

Many-body calculations similar to the calculations of energies described above show that 
q = 40500 MHz for Hg II and q = 956 MHz for Cs I. Note that A)( 133 Cs) = 2298157943 Hz. 
This is the exact value because the frequency 9192631770 Hz of the 6s F = 3 — 6s F = A hfs 
transition in 133 Cs, which is equal to AA (6s), is used as a definition of the metric second. 
Relative drifts in rates of atomic clocks based on Hg II and Cs will be 

Tt 4^4 = <Hg + ) - = (2.30 - 0.83)* = 1.47*. 

at A{LiS) Aq Aq a a a 



This is in good agreement with estimates based on the Fermi-Segre formula [ 2"2fl . The best 
sensitivity of various clock rate comparisons can be achieved when Hg + clocks are compared 
with an H maser (parameter q = for H). Assuming again at/at = 10 -15 yr _1 we will 
get 2.3 x 10~ 15 yr _1 for the frequency rate shift between H maser and Hg + clocks. Note that 
the ratio of hyperfine structure constants is also sensitive to the variation of nuclear magnetic 
(^-factors which may appear due to variation of strong interaction. 

One more interesting possibility is to use transitions between "accidentally" degenerate 
levels in the same atom. Such meta-stable levels exist, for example, in the Dy atom: two 
J = 10 opposite parity levels 4:f w 5d6s and 4f 9 5d 2 6s lying 19797.96 cm -1 above ground state. 
(This pair of levels was used to study parity non-conservation in Refs. [2"3[,2"5|]). There are 



other examples of "accidentally" degenerate levels in the rare-earth and actinide atoms and 
many close levels in other heavy atoms and ions (in the absence of degeneracy one should 
look for s — d or s — p transitions where the relativistic effects are larger). In the case of 
"accidental" degeneracy, the contribution of the relativistic correction to the frequency of 
the El transition in a heavy atom (~ 1000cm" 1 ) is compensated by the difference in the 
Coulomb interaction energies of the two configurations. However, if a varies with time, this 
compensation will eventually disappear. Thus, we have a correction ~ 1000 cm _1 ((^) 2 - 1) to 
the very small (< 0.01 cm -1 ) frequency of the transition. One can measure, for example, the 
time dependence of the ratio of frequencies for transitions between the hyperfine components 



of these two states. In the case of "accidentally" degenerate levels belonging to different 
electron terms in a molecule one can have enhanced effects of the change of both a and the 
nucleon mass. In the latter case the enhancement factor is the ratio of the vibration energy to 
the small frequency of the transition. The problem in these degenerate level cases is to find a 
transition with a small natural width. 
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TABLES 



TABLE I. Energy levels of C IV and Si IV with r 


espect to the continuous spectrum 


limit ( c m 1^ 

111111 U l L^lll j 


Ion 


State 


RHF 




Brueckner 


Experiment 


C IV 


2s x / 2 


519253 




520082 


520178 




2P1/2 


454054 




455508 


455694 




2p 3 /2 


453926 




455337 


455587 




3G?3/2 


195196 




195287 


195298 




3c?5/2 


195187 




195277 


195287 


Si IV 


3S1/2 


360613 




363840 






3P1/2 


290073 




292514 


292808 




3^3/2 


289606 




292036 






3G?3/2 


201807 




203480 


203721 




3d 5 /2 


201807 




203472 




a Reference 27 


|; numbers are 


rounded to the last digit before the decimal point. 




TABLE II. Ionization potential and excitation 


energies of Cr II (cm 1 ) 




State 


RHF a 


RHF + S b 


CI + E c 


CI + / x S d 


Experiment 






Ionization potential 






3d 5 6 S 5/2 


125889 


141067 


137208 


133815 


133060 




Excitation 


energies 






3d 4 Ap 6 F 3/2 


40230 


52943 


50233 


46777 


46905.52 


UHp 6 F 5/2 


40388 


53125 


50411 


46949 


47040.54 


3dHp 6 F 7/2 


40606 


53377 


50655 


47187 


47227.50 


3d 4 4p 6 P 3/2 


42422 


55171 


51902 


48466 


48399.19 


3d 4 4p 6 P 5/2 


42562 


55342 


52067 


48621 


48491.39 


3dHp 6 P 7/2 


42758 


55581 


52304 


48844 


48632.36 



a Single-configuration approximation. 

b As RHF but correlations with core electrons included in the second order. 

c Configuration interaction for 343 relativistic configurations; correlations with core are also included. 
d Screening parameter fd = 0.74 has been set for correlations of d-valence electrons with the core to 
fit ionization potential and excitation energies. 
Reference |7j. 



TABLE III. Dependence on a of the frequencies of the El atomic transitions of astronomic 



interest 


(cm x ). Here 


id = uj + qix + q 2 y where x = 


) 2 -i,y = 








Z 


Atom/Ion 


Ground state 


Upper states 




Qi 


<?2 


6 


C II 


2s 2 2p 


2 po 
M/2 


2s2p 2 


O/ Z 


74932. 617 a 


177 


3 








2s2p 2 


2 c 
^1/2 


96493. 742 a 


171 


3 










2s2p 2 


2 p 

1/2 


110625. l b 


173 


-3 










2s2p 2 


- r 3/2 


110666.3 b 


217 


3 


6 


C IV 


ls 2 2s 


2 Si/o 


\s 2 2p 


2 P , 
1 1/2 


64484. 094 a 


108 


8 








\s 2 2p 


r 3/2 


64591. 348 a 


231 


-8 


7 


N V 


2 s 


2 Si /o 
°l/2 




2 Pt /o 
^ 1/2 


80463. 211 a 


196 


-4 








2o 


r 3/2 


80721. 906 a 


488 


2 


8 


I 


2d 4 


3 P9 


2p z 3s 


3 Si 


76794. 977 a 


130 


-30 










2p z As 


3 Si 


96225.055 a 


140 


-20 


12 


Mg I 


3s 2 


So 


3s3p 


X P1 


35051.277 c 


106 


-10 


12 


Mg II 


3s 


2 Sl /9 


3p 


2 Pl /9 


35669.298 c 


120 











3p 


2p 


35760. 848 c 


211 





13 


Al III 


3s 


2 Sl /9 


3p 


2 Pl /9 


53682. 330 a 


216 











3p 


2 P 3 /2 

Of z 


53916. 598 a 


464 





14 


Si II 


3s 2 3p 


2 po 

M/2 


3s3p 2 


2 L>,/9 


55309. 352 a 


547 


-6 


14 


Si IV 


2p 6 3s 


^1/2 


2p 6 3p 


2p 

r 1/2 


71287.523 a 


362 


-8 








2p 6 3p 


2p 

Jr 3/2 


71748. 625 a 


766 


48 


20 


Cal 


As 2 




AsAp 


1 1 


23652.305 a 


300 





20 


Ca II 


3p Q As 


2q 

ill /o 


3v^Ad 


2p 

1/2 


25191. 512 a 


192 


16 










op <±p 


r 3/2 


25414.427 a 


420 


16 


24 


Cr II 


3d 5 


°5/2 


3d 4 Ap 


x 3/2 


46905. 17 d 


-1624 


-25 








3d 4 Ap 


r 5/2 


47040. 35 d 


-1493 


-21 










3J 4 4r) 


r 7/2 


47227. 24 d 


-1309 


-18 










3d 4 Ap 


6 Ps/2 

Of Z 


48398. 941 a 


-1267 


-9 










3d 4 Ap 


6p , 
X5/2 


48491. 105 a 


-1168 


-16 










3d 4 Ap 


6p , 


48632. 125 a 


-1030 


-13 


25 


Mn II 


3d 5 is 


7 s 3 


3d 5 Ap 


7 P 2 


38366. 184 a 


918 


34 










3d 5 Ap 


7 P 3 


38543.086 


1110 


19 










3d 5 Ap 


7 P 4 


38806.664 


1366 


Zl 


26 


be 11 


3d 6 As 


6 D 9/2 


3d 6 Ap 


6p.o 
*-^a /9 


38458.9871 


1449 


2 










3d 6 Ap 


6 

7/2 


38660. 0494 e 


1687 


-36 










3(fAp 


6rr 
r 11/2 


41968. 0642 e 


1580 


29 










3d 6 Ap 


6 F 9/2 


42114. 8329 e 


1730 


26 










3d 6 Ap 


6 F 7/2 


42237.0500 6 


1852 


26 










3d 6 Ap 


6 P 7 /2 


42658. 2404 e 


1325 


47 


30 


Zn II 


3d w As 


2 Sl/2 


3d w Ap 


2 Pl/2 


48480. 992 a 


1445 


66 








3d w Ap 


2 P 3 /2 


49355.027 a 


2291 


94 


32 


Ge II 


As 2 Ap 


2 P , 
M/2 


As 2 5s 


2 s V2 


62403. 027 a 


-575 


-16 



a Morton, Ref. g6 
b Moore, Ref. §7| 



c Pickering, Thorne and Webb, Ref. f28 |. 
d Sugar and Coliss, Ref. M. 



3 Nave et al, Ref. [30|. 



TABLE IV. Relativistic shift of energies of some methastable states of atoms which used as 
optical frequency standards (cm -1 ). 



Z Atom/Ion 


Ground state 


Upper states 




Qi 


Q2 


20 Ca I 


As 2 




AsAp 


3 Pi 


15210 


230 





38 Sr II 


5s 


2 Sl/2 


Ad 


2D 3/2 


14555.90 


2636 


96 






Ad 


2 D 5/2 






1 fin 


70 Yb II 


6s 


2 s 1/2 


5d 


2 D 3 /2 


24332.69 


9898 


1342 






5d 


2 D 5/2 


22960.80 


8298 


1570 


80 Hg II 


5d l0 6s 


2 Sl/2 


5d 9 6s 2 


2 D 5/2 


35514.0 


-36785 


-9943 






5d 9 6s 2 


D3/2 


50552.0 


-19377 


-12313 



